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We apply Mie scattering theory to study the interaction of magnetic spheres with microwaves in 
cavities beyond the magnetostatic and rotating wave approximations. We demonstrate that both 
strong and ultra-strong coupling can be realized for a stand alone magnetic spheres made from 
yttrium iron garnet (YIG), acting as an efficient microwave antenna. The eigenmodes of YIG spheres 
with radii of the order mm display distinct higher angular momentum character that has been 
observed in experiments. 
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I. INTRODUCTION 

Light-matter interaction in the strong coupling regime 
is an important subject in coherent quantum information 
transfer.Spin ensembles such as nitrogen-vacancy cen¬ 
ters may couple strongly to electromagnetic fields and 
have the advantage of both long coherence times^ and 
fast manipulation.® The “magnon” refers to the collec¬ 
tive excitation of spin systems. In paramagnetic spin en¬ 
sembles in an applied magnetic field, the spins precess 
coherently in the presence of microwave radiation, creat¬ 
ing hybridized states referred to as magnon-polaritons.®^® 

In the strong coupling regime coherent energy exchange 
exceeds the dissipative loss of both subsystems. The co¬ 
herent coupled systems is usually described by the Tavis- 
Cummings (TC) model,which defines a coupling con¬ 
stant g between the spin-ensemble and the electromag¬ 
netic radiation that scales with the square root of the 
number of spins. In ferro/ferrimagnets the net spin den¬ 
sity is exceptionally large and spontaneously ordered, 
which makes those materials very attractive for strong¬ 
coupling studies. The exchange coupling of spins in mag¬ 
netic materials also strongly modifies the excitation spec¬ 
trum into a spectrum or spin wave band structure. An 
ubiquitous experimental technique to study ferromag¬ 
netism is ferromagnetic resonance (FMR), i.e. the absorp¬ 
tion, transmission or reflection spectra of microwaves. In 
the weak coupling regime FMR gives direct access to the 
elementary excitation spectrum of ferromagnets,^^ includ¬ 
ing the standing spin waves in confined systems referred 
to spin wave resonance (SWR).^^ The strong coupling 
regime is studied less frequently, however, because the 
dissipative losses of the magnetization dynamics are usu¬ 
ally quite large. 

An exceptional magnetic material is the man-made 
yttrium iron garnets (YIG), a ferrimagnetic insulator. 
Gommercially produced high-quality spherical YIG sam¬ 
ples serve in magnetically tunable filters and resonators 
at microwave frequencies. By suitable doping becomes 
a versatile class of materials with low dissipation and 
unique microwave properties.^® YIG has spin density 
of 2 X 10^^cm“®,®^ and the Gilbert damping (recipro¬ 


cal quality) factor of the magnetization dynamics ranges 
from 10“® to 10“®,®®^®^ which facilitates strong cou¬ 
pling for smaller samples. Indeed, strongly coupled mi¬ 
crowave photons with magnons have been experimen¬ 
tally reported for either YIG films with broadband copla- 
nar waveguides (GPWs),®®^^'’ or YIG spheres in 3D mi¬ 
crowave cavities. ^®^^® A series of anticrossings were ob¬ 
served in thicker YIG films and split rings.The cou¬ 
pling of magnons in YIG spheres with a superconduct¬ 
ing qubit via a mircowave cavity mode in the quantum 
limit has been reported.^® An ultrahigh cooperativity 
C = g^ jKq > 10®, where k and 7 are the loss rates of the 
cavity and spin system, and multimode strong coupling 
were found at room^^ as well as the low^® temperatures. 

From a theoretical point of view, the standard TG 
model is too simple to describe the full range of coupling 
between magnets and microwaves. Also the rotating-wave 
approximation (RWA) (usually but not necessarily as¬ 
sumed in the TG model) is speaking applicable when 
the coupling ratio g/cOc ^ !> where Wc is the microwave 
cavity mode frequency. We may define different cou¬ 
pling regimes, ^^’^® viz. (i) strong coupling (SG) when 
0.01 < g/oJc ^ 0.1, (ii) ultrastrong coupling (USG)^® 
when g/ujc ^ 0.1, (iii) or even deep strong coupling 
(DSG) g/uJc ~ 1.^®" Gao et al.^ adapted the TG model to 
ferromagnets by formulating a first-principles scattering 
theory of the coupled cavity-ferromagnet system based 
on the Maxwell and the Landau-Lifshitz-Gilbert equa¬ 
tion including the exchange interaction. A effectively one¬ 
dimensional system of a thin film with in-plane magneti¬ 
zation in a planar cavity was solved exactly in the linear 
regime, exposing, for example, strong-coupling to stand¬ 
ing spin waves. Maksymov et alP'^ carried out a numerical 
study of the strong coupling regime in all-dielectric mag¬ 
netic multilayers that resonantly enhance the microwave 
magnetic field. A quantum theory of strong coupling for 
nanoscale magnetic spheres in microwave resonators has 
been developed in the macrospin approximation,^® but 
this regime has not yet been reached in experiments. 

Here we apply our classical method® to spherically sym¬ 
metric systems, i.e., a magnetic sphere in the center of a 
spherical cavity. This is basically again a one-dimensional 
problem that can be treated semi-analytically and has 
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other advantages as well, such as a homogeneous dipolar 
field and simple boundary conditions. The eigenmodes 
of magnetic spheres have been studied in the “magneto¬ 
static” approximation,in which the spins interact by 
the magnetic dipolar field, disregarding exchange as well 
as propagation effects, which may be done when A ^ a, 
where a is the radius of the sphere and A the wavelength 
of the incident radiation. Arias et treated the in¬ 
teraction of magnetic spheres with microwaves in the 
weak-coupling regime. In contrast, we address here the 
properties of the fully hybridized magnon-polaritons be¬ 
yond the magnetostatic approximation (but disregard the 
exchange interaction), including the propagation effects 
(reflection and transmission) of microwaves, thereby ex¬ 
tending the validity to A < a. We are admittedly still 
one step from the “exact” solution by disregarding the 
exchange (as treated and discussed by Cao et al.^). Our 
calculated microwave spectra are complex but help in un¬ 
derstanding some of the above-mentioned experiments. 

This manuscript is organized as follows. In Sec. II, we 
introduce the details of our model and derive the scat¬ 
tered intensity and efficiency factors for a strongly cou¬ 
pled system of a magnetic sphere and microwaves. In 
Sec. Ill, we present and discuss our numerical results 
that demonstrate the effects both due to the dielectric 
as well as magnetic effects on the scattering properties 
and compare our results with experiments. In Sec. IV, 
we conclude and summarize our findings. 


II. MODEL AND FORMALISM 

We model the coupling of the collective excita¬ 
tions of a magnetic sphere to microwaves in a spheri¬ 
cal cavity by the coupled Landau-Lifshitz-Gilbert and 
Maxwell equations. We employ Mie-type scattering the¬ 
ory, i.e., a rapidly converging expansion into spherical 
harmonics.We model the incoming radiation as 
plane electromagnetic waves with arbitrary polarization 
and wave vector that are scattered by a cavity loaded by a 
magnetic sphere with gyromagnetic permeability tensor 
In order to understand the experiments it is not 
necessary to precisely model the details of the resonant 
cavity. Instead, we propose a generic model cavity that 
is flexible enough to mimic any realistic situation by ad¬ 
justing the parameters. We consider a thin spherical shell 
of a material with high dielectric constant Cc/eo 3> 1, ra¬ 
dius R, and thickness <5 that confines standing microwave 
modes with adjustable interaction with the microwave 
source (see Fig. 1). The spherical symmetry simplifies the 
mathematical treatment, while the parameters R and <5 
allow us to freely tune the frequencies and broadenings 
of the cavity modes. 

The dynamics of the magnetization vector M is de¬ 
scribed by the LLG equation, 

= -yM X Heff + x (1) 


with a and 7 being the Gilbert damping constant and 
gyromagnetic ratio, respectively. The effective magnetic 
field Heff = He,t + a. comprises the external and 
(collinear) easy axis anisotropy fields Hoxt as well as the 
exchange field Hx = JV^M, with J being the exchange 
stiffness. Assuming that perturbing microwave magnetic 
field and magnetization precession angles are small: 

M(r, t) = Mg -I- m(r,t) (2) 

H(r,t)=Hext + h(r,t) (3) 

where is the saturated magnetization vector and m 
the small-amplitude magnetization driven by the rf mag¬ 
netic field h, we linearize the LLG equation to 

dtm = -jMs X - ^dtm^ - ym x (4) 

where = H^xt and -|- h. The response of 

ferromagnetic spheres is affected by exchange when their 
radii approach the exchange length. Since the latter is 
typically a few nm, we hereafter disregard the exchange 
interaction and concentrate on the dipolar spin waves. In 
the frequency domain and taking the z direction as the 
equilibrium direction for the magnetization, 

iujm = z X (clImL — wnm -I- icoam), (5) 

with wm = yAfg and wh = iHo- We may recast Eq. ( 6 ) 
into the form m = ^ • h.The magnetic susceptibility 
tensor ^ is related to the magnetic permeability tensor 
by V = Mo(^ + V)- We find 

/l + X 0\ 

V = Mo 1 + X 0 (6) 

Vo 0 1/ 

where y and k are given by. 


(wh — iauj)uJM 
^ (wh — 

(7) 

OJUJm 

(8) 

(wh — iOLLoY — 


The permeability tensor appears in the Maxwell equa¬ 
tions for the propagation of the electromagnetic wave in 
a magnetic medium. 

Inside a spatially homogeneous medium a monochro¬ 
matic wave with frequency w, 

V X E = iwb, V X h = —iwD (9) 

V • D = 0, V - b = 0. (10) 

The constitutive relation between the magnetic induction 
b, electric displacement D, magnetic field h, and the elec¬ 
tric field E inside this medium are 

b = V • h, D = CspE. 


( 11 ) 
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where egp is the scalar permittivity of the medium. It fol¬ 
lows from Eqs. (10) and (11) that the magnetic induction 
b satisfies the wave equation, 

V X V X (mo ■ b) - /c^pb = 0 (12) 

with = w^espMo- 

The surrounding (nonmagnetic) medium is homoge¬ 
neous and isotropic with scalar magnetic permeability mo j 
divergenceless magnetic field, and simplified wave equa¬ 
tion V^b + fcgpb = 0. Due to the spherical symmetry it is 
advantageous to expand the magnetic field h into vector 
spherical harmonics 


h = ^?7r, 


r) -f ] 


(13) 


where n runs from 1 to oo, and m = —n, ■ ■ ■ ,n with 
prefactors fj^m = ??nmfco/(wMo), 


Vnm — ^ -^0 


2n-|-l {n — m)l 


1 1/2 


n{n -I- 1) (n -I- m)\ 


(14) 


Eq is the electric field amplitude of the incident wave . 
The vector spherical harmonics 

MW(fc,r) = zW(A:r)X„^(f), 

fcNW(fc,r) = VxM«(fc,r). (15) 


Zn\kr) are spherical Bessel functions, X„m(f) = 
LEnm(r)/\/i2(i2 + 1) Spherical harmonics and L = —ir x 
Vr the angular momentum operator with Vr the gradi¬ 
ent operator. The electric field distribution is obtained by 
E = (i/wc)V X h. By invoking the vector spherical wave 
function expansion for b and • b in the wave equa¬ 

tion Eq. (12) leads to the dispersion relation for k (w). 

We match the field distributions inside and outside the 
cavity to obtain the scattering solution for incident plane 
microwaves. The field inside the spherical shell must be 
regular, while the scattered component has to satisfy the 
scattering wave boundary conditions at infinity. These 
conditions are fulfilled by adopting the first kind of spher¬ 
ical Bessel function jn{x) as radial part for the internal 
distribution and the first kind of spherical Hankel func¬ 
tion hn'^ (x) for the scattered component outside the cav¬ 
ity 

bs = ^ ( Vnm 
nm 

(16) 

The unknown scattering coefficients Cnm and dnm are de¬ 
termined by the boundary conditions at the interface. We 
consider here the situation in which the magnetic sphere 
is illuminated by a plane wave with arbitrary direction of 
propagation and polarization as indicated in Fig. 1. The 
incident field can be expanded as 


CnmX^2i(^0i r) -|- dnm'^^nmikoj r) 


^ ^ Vnm 


^nm NW(fco,r) + '^nm'^^hmiko^ I*) 


( 17 ) 



Figure 1. (Color online) Plane wave with wave vector ko com¬ 
ing in at an arbitrary angle hits a large spherical cavity mod¬ 
eled by a dielectric spherical shell of radius R, thickness 5, and 
permittivity tc- The spherical cavity is loaded with a magnetic 
sphere of radius a centered at the origin of the coordinate sys¬ 
tem. 


The expansion coefficients Umn and Vmnj 

Unm = [per„m(cos0fe) - ip 4 ,^nm{cos 0 k)] , (18) 

Vnm = [PeTinm{cOs9k) - /p0f„m(cOS 6»fc)] (19) 

contain all information about the polarization vector and 
direction of propagation, where p = (peOk +P 44 >k) is the 
normalized complex polarization vector, with |p| = 1 and 
9k{4’k) is the polar (azimuthal) angle of kg. Two auxiliary 
functions are defined by 

Unm ~ inm^. (COS^), Tnm ~ inm^^Z^n (ct)S^), 

Sin a da 

( 20 ) 

with tnm = i '^Vnm/Eo and P™{x) the first kind associ¬ 
ated Legendre function. 

In order to solve the full scattering problem including 
the cavity we match the fields outside the cavity caused 
by the incoming plane microwave and the spacer region 
separating the magnetic particle and cavity. In the latter, 
spherical Bessel functions of both the first and second 
kind have to be included into the expansion. At the sur¬ 
face of the magnetic sphere (r = a) we adopt the stan¬ 


dard boundary conditions 

hj X Gj. = hmid ^ (^^) 

Ej X Gj. = Einiid ^ (^2) 

while at the surface of the cavity, assuming that its thick¬ 
ness is much smaller than the wavelength, 

[bmid X Gj. = ^ [Gj. X Eo^i] X Gj., (2^) 

Eimid X Gj. = E/Qnt X Gj.. (^4) 


nm 
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Figure 2. (Color online) Scattering intensity |Sip as function 
of scattering angle 0 and frequency oj/2ti is shown for (a) a 
dielectric sphere of radius a = 1.25 mm and relative permit¬ 
tivity e/eo = 15 and for (b) the same sphere in a cavity of 
radius R — 1.6 mm. In (c) the scattering intensity is plotted 
for the same cavity as function of frequency and loading rate 
a/R. The dashed lines are guides for the eye. 



H„/M 
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Figure 3. (Color online) Panel (a) shows the scattering effi¬ 
ciency factor Qsca as function of normalized magnetic field 
Ho/Ms and frequency aj/27r for a YIG sphere of radius a = 
2 mm and relative permittivity e/eo = 15. Panel (b) shows 
results for a non-magnetic dielectric sphere. The character of 
the microwave modes sufficiently far from the anti-crossing 
with the spin waves is labeled by the spherical harmonic in¬ 
dices (n, m). 


The indexes mid and out indicate the regions within and 
outside of the cavity, respectively. The unit vector is 
the outward normal to the surfaces and ^ = zw(ec — eo)<5 
with permittivity of the cavity shell e^- By matching the 
field distributions in the different regions the scattering 
coefficients are determined, from which we calculate the 
observables. At distances sufficiently far from the cavity, 
i.e., in the far field zone, the intensity of the two polar¬ 
ization components Ii and I2 are 


h 

h 




(25) 

(26) 


where 0 {(/>) is the polar (azimuthal) angle of the observer 
at distance r. The scattering amplitude functions are 


Si{9,(j)) = ^ [dnmfnm.{cOs9) + Cnm^nm{cOs9)] 

nm 

(27) 

S2{9, (t>) = ^ [dnm^nm{.COs9) + CnmTnm{cOS 9)] 

nm 

(28) 


where the coefficients Cnm and dnm characterize the scat¬ 
tered component of the fields outside the cavity. We may 
now compute the scattering and extinction cross sections 
as well as their (dimensionless) efficiencies Qsca and Qext, 
which are the cross sections normalized by ttR^, the geo¬ 
metrical cross section of the cavity: 


Qsca — R/^ ^ ^ (|Onm| T |^nm| ) j 

^ nm 

Qext = ^ 2^2 ^ + "^nm^nm) ■ 


(29) 

(30) 


The extinction cross section represents the ratio of (angle- 
integrated) emitted to incident intensity, i.e., with and 
without the scattering cavity/particle between source 
and detector. This factor measures the energy loss of the 
incident beam by absorption and scattering. The series 
expansion in Eqs. (27)-(30) is uniformly convergent and 
can be truncated at some point in numerical calculations 
depending on the desired accuracy. In the next section we 
present our results with emphasis on the dielectric and 
magnetic contributions to the microwave scattering. 


III. RESULTS 

Here we present numerical results on the coupling of 
microwaves with a ferro- or ferrimagnet in a cavity based 
on our treatment of Mie scattering of the electromagnetic 
waves as exposed in the preceding section. It applies to 
a dielectric/magnetic sphere centered in a (larger) spher¬ 
ical cavity, but both may be of arbitrary diameter oth¬ 
erwise. We are mainly interested in the coherent cou¬ 
pling between the magnons and microwave photons in 
the strong or even ultrastrong coupling regimes that can 
be achieved by generating spectrally sharp cavity modes, 
by increasing the filling factor of the cavity, or simply 
by increasing the size of the sphere. The RWA, however, 
tends to break down as the coupling increases. This has 
led to different coupling regimes beyond the weak cou¬ 
pling, TC region, i.e., strong (SC) and ultrastrong (USC) 
coupling regimes. In the SC region coupling strength has 
to be comparable or larger than all decoherence rates, 
while in the USC it has to be comparable or larger than 
appreciable fractions of the mode frequency, g/coc ^ O-l- 

We adopt the forward scattered intensities Ii ^ 
|S'i(6* = 7r/2,(/) = 7r)p and scattering efficiency fac¬ 
tors as convenient and observable measures of the mi¬ 
crowave scattering by a spherical target. In order to com- 
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Figure 4. (Color online) Scattering efficiency factor Qsca plot¬ 
ted as function of normalized magnetic field Hq/Ms and fre¬ 
quency uj/2tt for a YIG sphere of radius a = 1.25 mm and 
relative permittivity e/eo = 15 (a) in the center of a spherical 
cavity of radius R = 1.6 mm and (b) without cavity. 


Figure 5. (Color online) Scattering efficiency factor Qsca as 
function of normalized magnetic field Hq/Ms and frequency 
tj/27r for a YIG sphere of radius a = 1.25 mm and relative 
permittivity e/eo = 15 (a) in the center of a spherical cavity 
of radius R = 1.6 mm and (b) in the absence of the cavity. 
Dashed lines indicate the frequency range in Figure 4. 


pare results with recent experiments, we chose param¬ 
eters for YIG with gyromagnetic ratio 7/(27r) = 28 
GHz/T, saturation magnetization^^ ITqMs = 175 mT, 
Gilbert damping constant^®^^^ a = 3 x 10“^, and rela¬ 
tive permittivity"^^ e/co = 15. Without loss of general¬ 
ity we consider microwaves incident from the positive 
X direction (9k = 7r/2 and (pk = 0) and polarization 
(pg,P(p) = (1,0), so its electric/magnetic components are 
in the —ziy directions (static magnetic field and magneti¬ 
zation Ho||z). Forward scattering is monitored by setting 
9 = -K12 and (j> = tt in Eq. (27). We also explore the de¬ 
pendence of the observables on the scattering angles. We 
can remove the cavity simply by setting ^ = 0. 

In Fig. 2 the scattered intensity |5'i(0,7r)p is depicted 
as a function of frequency a;/27r and scattering angle 
9 focusing first on a non-magnetic sphere with radius 
a = 1.25 mm. The angular dependence of the scattering 
with and without a cavity (with R = 1.6 mm) is plot¬ 
ted in panel (a) and (b), respectively. The eigenmodes 
of the dielectric sphere show s, p and d-wave charac¬ 
ters in Fig. 2(a). s-wave scattering dominates as long as 
the wavelength (reduced by Csp) does not fit twice into 
the sphere, i.e., A ^ a^/esp/eo- The spherical cavity, on 
the other hand, limits the isotropic scattering regime to 

A ~ R\/ ^sp/Co- 

In Fig. 2(c) we plot the forward scattered intensities Ii 
as function of the load of the cavity by a dielectric sphere. 
The eigenfrequencies of the cavity remain constant, while 
those confined to the sphere shift to lower frequencies as 
~ a~^. At high loading rate the cavity modes are strongly 
mixed with the modes in the sphere and all of them bend 
towards lower frequencies. 

Magnetism of the spheres can affect the microwave 
scattering properties strongly, but the issue of hybridiza¬ 
tion of cavity and sphere resonant microwave modes is 
still present. A sufficiently large YIG sphere alone can 
therefore provide strong coupling conditions to the mag¬ 
netization even without an external resonator. To this 
end, the linear dimension of the YIG sphere must be of 
a size that allows the internal resonances of the sphere 


to come into play in the microwave frequency range, i.e. 
when ka ^ Tr-^/eo/esp or A ^ 2ay/Csp/eo- We therefore 
have a (narrow) regime a-^/csp/eo ^ ^ ^ 2a-\/esp/eo or 
7.75 mm ^ A ^ 15.49 mm (for Fig. 3) in which strong 
coupling and s-wave scattering can be realized simulta¬ 
neously without a cavity. YIG spheres can typically be 
fabricated with high precision for radii in the range^^ 
a = 0.9 — 2.5 mm. In Fig 3 for a a = 2 mm YIG 
sphere we observe a strong anticrossing between the lin¬ 
ear spin wave modes and the sphere-confined standing 
microwaves. The YIG sphere is therefore an efficient mi¬ 
crowave antenna that achieves strong and ultra-strong 
coupling without a cavity. It should be noted that previ¬ 
ous works,which have revealed the possibility to use 
all-dielectric as well as all-magneto-dielectric resonators 
without external resonator, were not considered strong 
coupling. 

Our results help to interpret recent experimental re¬ 
sults on YIG spheres in microwave cavities with reported 
coupling strength that are comparable with the magnon 
frequency,i.e., in the ultrastrong-coupling regime. In 
Fig. 4 the scattering efficiency factor is shown as a func¬ 
tion of Hq/Ms and a;/27r. Panel (a) addresses a YIG 
sphere of radii a = 1.25 mm in a spherical microwave cav¬ 
ity of radii R = 1.6 mm, chosen to be close to the lead¬ 
ing dimensions of the cavity in the experiments. Panel 
(b) holds for the same YIG sphere but without cavity. 
The obvious anticrossing in Fig. 4(a) is a signature of 
the emergence of the hybrid excitation that we refer 
to as magnon-polariton. The anti-crossing modes are la¬ 
beled by the mode numbers (n,m). For given n there 
are two m = ±n anti-crossing modes with coupling 
strengths > gn,-n, where gn,m is the effective cou¬ 
pling strength of the magnon mode (n, m) to the cavity. 
Fig. 4(a) indicates that the ultrastrong-coupling strength 
is indeed approached since a splitting of g/2Tr = 2.5 GHz 
is achieved at a resonance frequency of a;/27r ~ 37.5 GHz. 
Beside the main anticrossing with the (2, 2) and (2, —2) 
cavity modes, we observe tails from other anticrossings 
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with the (3,3) and (3,-3) modes at higher frequencies, 
as well as the (1,1) and (1,-1) modes at lower frequen¬ 
cies, which are standing electromagnetic resonance modes 
confined by the YIG sphere. We may interpret these as 
nearly pure spin wave modes that acquire some oscilla¬ 
tor strengths by mixing from far away resonances due to 
the ultrastrong coupling with standing microwaves. This 
can be verified by checking the scattering efficiency fac¬ 
tor in the absence of the cavity as in Fig. 4(b), which 
emphasizes the antenna action of the YIG sphere. 

Zhang et al.^^ indeed report additional, weakly- 
coupled “higher modes,”but without explaining their na¬ 
ture. They report ultrastrong-coupling between magnons 
and the cavity photons only in the frequency range of 
35 — 40 GHz, but data at lower frequencies are not given. 
In Fig. 5 we extend the plots in Fig. 4 to a larger fre¬ 
quency interval. We observe that the main anticrossing 
in the frequency range of 35 — 40 GHz is caused by the 
n = 2 modes, while hybridized modes originating from 
the n = 1 resonance exist at the lower frequencies. The 
unperturbed modes between the anticrossing gaps are 
therefore not only due to the higher modes, but lower 
modes with n = 1 also contribute by the ultrastrong¬ 
coupling. Two significant curves in the left and right side 
of the higher unperturbed modes originate from the an¬ 
ticrossing modes n = 1 (the left one) and n = 2 (the 
right one) of the YIG sphere itself, as is more clear in 
Fig. 5(b) (the computed lines are broader because we 
use a relatively large k for computational convenience). 
We thereby find again that the strong-coupling magnon- 
polariton may form also without cavity. 

We concentrated on the dipolar spin wave excitations 
driven by magnetic fields that are strongly inhomoge¬ 
neous due to a large dielectric constant. We disregard 
here exchange interactions, thereby limiting the validity 
of the treatment to YIG spheres much larger than the so- 
called exchange length that for YIG is only a few nanome¬ 
ters. In other words, we cannot properly describe all spin 
waves with relatively large wave number or frequencies 
relatively much higher relative to the FMR frequency. In¬ 
deed, in the planar configuration spin wave resonances 
are observable for rather thick films.® Exchange-induced 
whispering gallery modes on the surface of the YIG might 
therefore be observable even in thicker spheres, but their 
treatment is tedious and beyond the scope of the present 
paper. 


IV. CONCLUSION 

In this paper we implement Mie scattering theory to 
study the interaction of dielectric as well as magnetic 
spheres with microwaves in cavities by the coupled LLG 
and Maxwell equations, disregarding only the exchange 
interaction. We are mainly interested in the coherent 
coupling between the magnons and microwave cavity 
modes in the strong- or even ultrastrong-coupling regimes 
characterized by the mode-dependent coupling strengths 
gn,m- We reveal that while in the presence of a spheri¬ 
cal cavity both strong and ultrastrong coupling can be 
realized by tuning the cavity modes and by increasing 
the filling factor of the cavity. Surprisingly, these regimes 
can also be achieved by removing the external resonator, 
due to the strong confinement of electromagnetic waves 
in sufficiently large YIG spheres. In this regime, higher 
angular momentum eigenmodes of the dielectric sphere 
participate and the scattering shows s- as well as p-wave 
character. We thereby transcend studies that focus on 
dipolar spin waves in a magnetostatic framework®^’®^ 
by considering propagation effects via the full Maxwell 
equation. Our study might be useful in designing opti¬ 
mal conditions to design cavities in which YIG spheres 
are coherently coupled to, e.g., superconducting qubits, 
in microwave cavities for coherent quantum information 
transfer.^® 
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